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1. Introduction 

Studied for the first time in 1925 pQ, the Ising model is one of the most important models 
of statistical physics. The two dimensional case, solved exactly for the first time in 1944 
by Onsager j2j, is the prototype of systems with second order phase transition and 
non Gaussian critical exponents. It has therefore been studied extensively by various 
exact and approximate methods. One important result is when Schultz, Mattis and Lieb 
[31 IU E] used transfer matrix method, Jordan- Wigner transformation and fermionization 
to simplify the Onsager solution in zero magnetic field, whereas fruitful links with ID 
Quantum Field Theory and Conformal Field Theory have been developed in a more 
recent period to provide a more general frame for studying 2D critical systems. In 
particular, a way to express the Ising Hamiltonian as a Gaussian Grassmannian action 
was established long time ago [HI Ej , and this idea was extended by Plechko to compute 
the partition function of the 2D Ising model in zero field for a large class of lattices 
[HUH], using operator symmetries that simplify the algebra of transfer matrix. 
The 2D Ising model has also been a starting point to study some boundary problems, 
for example the effect of a boundary magnetic field (BMF) on the propagation 
of a domain wall which is similar to a wetting or pinning problems. It is quite 
remarkable that exact results can be obtained in this case since the model with 
a general uniform field has not yet been solved except at the critical temperature 
[TUJ [TT] where conserved quantities have been found. In a series of papers, McCoy 
and Wu computed the partition function associated with a uniform BMF on one 
side of a square lattice and evaluated the boundary correlation functions, using dimer 
statistics [12 H3J EH E3 EE El With transfer matrix and spinor methods, some exact 
results have also been found for configurations with different fixed boundary spins or 
equivalently infinite BMF |T%1 IT9"1 !2Uj. or finite BMF |2*T] . These methods use mainly 
Jordan- Wigner transformation and spinor method related to rotation matrices [221 l2Hj- 
One interesting case is a configuration with infinite BMF on the two sides of an infinite 
strip with opposite sign|18j. A single wall domain pinned along the middle of the 
strip is present at low temperature, separating two regions of opposite magnetisation. 
Diffusive or sharp interfaces occur depending on temperature. Another case is when 
each boundary is composed of spins up followed at one point by spins down alongside 
of an infinite strip [2U [23 • From this point, an interface develops perpendicular to 
the strip inside the bulk. Other possible configurations with fixed boundary spins were 
studied using the same methods [201 [2H] • A solution with finite BMF was also proposed 
using Boundary Quantum Field Theory [22j, or Conformal Field Theory in continuous 
systems [28 J. 

Our aim here is to present an alternative derivation of the exact partition function 
including a general BMF, extending the method introduced by Plechko. This method 
appears to be simpler than the McCoy and Wu's derivation, or transfer matrix methods, 
and allows exact calculations for more complicated cases of BMF with finite amplitude. 
By this we demonstrate that the partition function can be simply expressed as a 
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Grassmann path integral of a ID Gaussian action with general or random magnetic fields 
after integrating over the bulk degrees of freedom. In the special case of an homogeneous 
BMF, we are able to compute the free energy on the lattice, and boundary spin-spin 
correlation functions. We check that our formula is equivalent to the one obtained by 
McCoy and Wu [121 HE]. Moreover, we also apply this method to obtain the exact 
partition function and the free energy when there is an interface developing between 
two finite opposite fields h and —h on one boundary side, as the simplest application 
of an inhomogeneous BMF. To our knowledge, this case has not been considered in 
the literature, except for infinite fields [TH] , and is very similar to a problem of an 
interface pinning on a boundary with the strength of the pinning that can be tuned 
(here the magnitude of the field). This method can be easily extended for more general 
configurations and therefore is useful for instance to studying wetting problems [T8|l29j. 
The article is organised as follow: In section 2, we introduce the notations used 
throughout the article. In section 3 and 4, following Plechko's method, we obtain a 
Grassmann path integral representation of the partition function and the main action. 
In section 5, as a preliminary and useful exercise, we solve in the same way the ID Ising 
model with an homogeneous magnetic field, in order to introduce the method to the 
reader. Section 6 is dedicated to the explicit calculation of the partition function for 
the two dimensional Ising model; we then give the corresponding boundary ID action 
after integrating over the bulk variables. At this stage, we compare our results with 
those of McCoy and Wu, and Au-Yang and Fisher, taking the thermodynamic limit. In 
the section 7, we compute the two-point correlation function on the boundary, which is 
necessary, in section 8, to obtain the expressions of the partition function and the free 
energy in the case of a pinned interface on one boundary with an inhomogeneous BMF. 



2. General notation 



In the following we consider the Ising model on a square lattice of size L with spins 
a mn — ±1. For sake of simplicity, we limit ourself here to the case where the coupling 
constant J is the same in both directions. The method works however if there are two 
different coupling constants along vertical and horizontal bounds. Until section 4 we 
consider inhomogeneous magnetic fields h n , placed on the sites of the first column m — 1 
(see Figure [TJ. Periodic boundary conditions for the spins are imposed in the direction 
parallel to the magnetic field line, a m \ = & m L+i, and free boundary conditions in the 
transverse direction, formally equivalent to a^ n = (TL+in = 0. 
The Hamiltonian is then given by 

L L 

T~t — — J } {Omn&m+ln + &mn&mn+l) ~ ^ h n Oi n . (1) 
m,n=l n=l 

The partition function (PF) Z is defined as 
Z = exp(-/3W), 



Figure 1. Description of the model on a periodical lattice with free conditions in one 
direction. 



where (3 = 1/ksT and the sum is over all possible spin configurations. We can write, 



using a r 



[cosh(/5J)] 2i Y\ cosh(^/i n )(l + uain) ■ JJ (1 + ta mn a m+ln ){l + to 

mn&mn+l ) 



n=l m,n=l 

with u n = tanh(/3/i n ) and t = tanh(/3J). We then define Q as 



Q[h] 



Tr 

O'mn 



J^J (1 + ta mn a m+ i n )(l + ta mn a mn+ i) ■ + u n a r . 



m,n=l 



n=l 



(2) 



where Tr^ is the normalised sum |^ cr= _ tl - This PF has already been calculated by 
McCoy and Wu ^2] for a uniform field on a boundary. This was done in two steps. 
First they proved that this PF is the Pfaffian of a matrix, using dimer statistics, they 
then performed the direct calculation of this Pfaffian. Here we treat the problem in 
a different way. Our idea is to generalise the elegant method introduced by Plechko 
[HI E] for the 2D Ising model in zero field since this method appears to be simple, and 
provides a direct link with Quantum Field Theory. It is indeed straightforward to obtain 
the expression of the quadratic fermionic action. Here we show that we can derive such 
a quadratic action in presence of a general boundary magnetic field. 



3. Transformation of the PF using a Grassmann representation 

Following Plechko, we introduce pairs of Grassmann variables in order to remove the 

local interaction between spins. We briefly define some useful tools using Grassmann 

algebra. For more details, we refer the reader to the book by Nakahara [3T)j . 

A Grassmann algebra A of size N is a set of iV anti-commuting objects {a,}^^ 

satisfying: 

V 1 < i,j < N, didj = —ajdi, 



which implies af = 0. Functions defined on such an algebra are particularly simple, they 
are polynomials. It is possible to define a notion of integration [30J with the following 
rules: 



da a = 1, / da 1 = 0, (3) 
and for any function f(a), 

da Ha) = ?M. (4) 
With these definitions, Gaussian integrals are expressed by 

/ Y[ da*dai exp ( a^a* ] = det M. (5) 

We also define a trace operator over an algebra A = {a, a*} as 

%t If (a, a*)} = [ da* da /(a, a*)e aa *, (6) 

a,a* J 

with the simple rules 

Tr [1] = 1, It [aa*] = 1. 

a, a* a, a* 

This operator will be useful in the following, and the subscripts may be omitted 
implicitely when the trace is performed over the Grassmann variables that are present 
in the expressions inside the brackets. Grassmann variables are introduced in the PF 
in order to decouple the spin variables. Terms containing the same spin are then put 
together and the sum over the spin configurations is performed. We use the fact that 

1 + too' = J da*da(l + aa) (1 + ta*o')e aa * = Tr[(l + aa) (1 + to* a')} (7) 

and follow closely the notation in reference jH] . We consider the following link variables 

V'mn = 1 + iCmn^m+ln) , . 

/(2) 1,4 W 

In order to decouple the products of two spins, we can express each object as a trace 
operator over a product of two Grassmann polynomials using equation (|7j) 

^(2) -Tr\B B* ] 

where 

■Amn 1 ~t~ (^mn^mni 1 -|- ta m _^ n O mn , (10) 

B mn 1 "I - b mn a mn} B mn 1 -|- tb mn _^o mn . 

The next step is to gather the different terms corresponding to the same spin. For 
that we consider the mirror symmetries introduced by Plechko and the operations 
of moving Grassmannian objects. Such operations are possible only within the trace 
operator. For example, the functions ipmn and ipmn are real functions, and therefore 
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commute with each other. We can also write 4>mn^mn = Xt [(A mn A! m+ln )(B mn B^ nn+1 )\ . 
Inside the brackets [...], the two groups (A mn A* m+ln ) and (.B m „.B^ n+1 ) can be moved 
independently since we can perform the trace on each of them separately, which gives 
real and therefore commuting quantities like ifiml, or tprnn that can be moved and inserted 
everywhere, and, at the end of these operations, we reintroduce the different integrations 
or trace function. For example we have 

^ (1 U (2) (11) 
= Tr [(A Tnn A m+ln )(B Tnn B rnn+1 )~j = Tr [(B mn B mn+1 )(A mn A m+1 
= Tr [(B mn (A mn A m+ln ) B mn+1 )~j = Tr ^( y A mn (B mn B mn+1 ) A m+ln j 

where the brackets (. . .) define the commuting group of Grassmannian objects. In 
general, if we have three groups of commuting objects {(OiO*)} i= i,, 3 we easily obtain 
the following mirror symmetries that can be applied to the objects (110)1 : 

{OlO^OlO^OtOz) = (^(0 2 *(0*0 3 )0 2 )Oi) (12) 

It is also important to treat the spin boundary conditions separately from the bulk 
quantities in order to obtain an expression valid not only in the thermodynamic limit, 
but for all finite values of N. Indeed, in the direction parallel to the line of magnetic 
fields where <J m L+i = cr m i, the corresponding link element can be expressed as 

^l = St [B mL B* mL+1 ] = Zt [B* ml B mL ] , (13) 

where B* ml = 1 + tb* mG a m i. The equality (fTB^) associated with definitions (fTUj) imposes 



—b* mL . The periodic conditions on spins therefore lead to anti-periodic conditions 



h* 

on Grassmann variables. In the transverse direction, we have o"o n = <JL+in — 0, 
corresponding to free boundary conditions. This implies the boundary conditions on 
Grassmann variables a^ n = and therefore A* ln = 1. The PF Q[h] (J2J) can be written 
in terms of the ipmn as 

L L 



Q[h] = Tr 



m,n=l 



n=l 



(14) 



Using the mirror symmetry ()12|). the boundary terms can be written as 



n=l 

n* 



(i) 

Ln 



(2) 
mL 



n-4 

n=l 
L 

J\ B ml 



m=l 



m=l 



IP 

m=l 



mL 



(15) 



where we introduce the notation 



N 



k=l 



A 1 -A 2 ---A h 
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N 



Y[A k = A N -A N _ 1 ---A 1 . 



k=l 



These products can be reorganised as follow 0: 



1 «-^ = s 

n=l m=l 



m=l n=l m=l 

For the bulk elements, we obtain the following arrangement 



U B *ml H A Ln X\B mL 



(16) 



m=l 

L-l L 



mn+1 



m=l 



m=l 



n=l m=l m=l 

nn«-n^ (i) 



n^i n nve n s 

m=l n=l m=l m=l 

L* / L L 



L-l L 



J Ln 



It 



n=l m=l 



71=1 



n=l \ m=l 



m=l 



(17) 
(18) 
(19) 



where we use the fact that ^in are commuting objects as well as the product 
reorganisation: 



o* 1 (o 1 o* 2 )(o 2 o* 3 )o 3 = l[o* m o r , 



(20) 



m=l 



We now insert the product over the remaining ipm n inside the previous expression (jTHj) 



L-l 



m=l 



m=l 



Bl n A ln ■ Y\_A* nn B* nn A Ti 



m=2 



(21) 



and finally obtain 



n ^(1)^(2) = 



m,n=l 



H^ln^l" I Y^AmnBmnAmn 4 Y\_ B-mn J B\ n 



71=1 



, m=2 



m=2 



(22) 



The PF is rewritten in a way that the sum on each spin can be performed by iteration 



Q[h] = Tr Tr 



l[B* ln A ln (l + u n a ln ) l[A 

mnBrnnAmn ' J^J B mn I B\ n 



71=1 



, m=2 



m=2 



(23) 



In fact, as we do not yet fermionize (1 + M n cr ln ), we reproduce here the Plechko's 
derivation in the special case of free-periodic boundary conditions 0. This expression 
is the basis of the rest of this paper. The sum over the spins outside the magnetic 
field region will lead to a quadratic action over Grassmann variables which therefore 
commute with the rest of the elements belonging to the first column (l,n). 
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4. Grassmannian representation of the action on the lattice 

The key point in equation is that the trace over the spin configurations is performed 
in an iterative way. For example, the first summation is done on spins a^n in the products 
A* Ln B* Ln AL n B in- This operation leads to a quantity which is quadratic in Grassmann 
variables, which can be put outside the general product (J2BJ). The same operation is 
then performed on spins er^-in and so on. This makes Plechko's method efficient for 
the 2D Ising model in zero field. With a uniform magnetic field in the bulk, the spin 
trace over the product of the four previous operators would lead to a quantity which 
is linear and quadratic in Grassmann variables and does not commute with the other 
products. However a BMF affects only the last products depending on spins a% n and 
this makes the problem very similar to a ID Ising model in a uniform magnetic field: 
It is therefore solvable. 



4-1. Trace over spins inside the bulk 

For spins cr mn inside the bulk, 1 < n < L and 2 < m < L, we have to evaluate 
successively 



with 



T r [^mn^rnn^-mnBmn] — exp(Q mn ), 



a mnbmn + t a m-ln^mn-l + ^( a m-ln ^mn-l)( a mn + b mn ). 



(24) 



These terms commute with all Grassmannian terms and can be pulled out of the 
remaining products. We obtain 



Tr 

\S'mn\m—1..L 



m=2 



m=2 




(25) 



The PF can now be written as 



Q[h] = Tr 



exp 




FT Tr ((1 + u n a ln )B* ln A ln B ln ) 

m „ 



01 

n=l s_ 



boundary spins d\ n 



(26) 



4-2. Trace over the boundary spins 

In the expression (j2*o]) . we can evaluate separately the trace over the spins <Ti n , leading 
to 



Tr[(l + u n cr ln )Bl n Ai n Bi n } = 1 + ax n b\ n + t&tn-i( a in + hn) + u n ^n + Wn*&* n -iain&in, 

h n = a ln + b ln + tb* ln _ v (27) 

The presence of a magnetic field on a site introduces a linear Grassmann term. This 
term no longer commutes with the others, and we need to compute the product in (|2*o]) 
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carefully. 

We would like to change artificially the fixed boundary conditions to periodic ones in 
order to simplify the subsequent calculations based on Fourier transformation. The 
quadratic part of equation (|2*7j) is equal to Qi n with fixed boundary conditions: We can 
write 



hn = Qi n -tal n h n , 



(28) 



Qln = ainhn + t^lnKn-l + + b ln~l)( a ln + &ln), 

where we introduce the boundary quantities a* 0n = a* Ln . Q\ n corresponds to periodic 
boundary conditions for the Grassmann variables (or anti-periodic conditions for the 
spins). We obtain 

Tr [(1 + u n a ln )B* ln A ln B ln ] = exp (Q? n + u n h n - ta* Ln h n ) . (29) 

flu 

The correction to periodic conditions due to the free boundary conditions for the spins 
is included in — ta* Ln L n . 



4-3. Grassmann variables associated with the magnetic field 

Here we introduce a pair of Grassmann variables (h n , ft,*) associated with the BMF. In 
the rest of the article, we will refer to it as the fermionic magnetic field. We have 



exp(« n L n ) = 1 + w n L r , 
therefore 



dh* n dh n (1 + u n h n ){l + h* n h n )e 



llnht 



(30) 



Tr [(1 + u n a ln )B* ln A ln B ln ] = Tr exp Q p ln + (h* n - ta* Ln )h n + u n h n 
We now perform a translation in the fermionic magnetic field 



TT I, TJ* V* +„* 



(31) 



which leads to 

L 



TT Tr [(1 + u n a ln )B* ln A ln B ln ] = Tr 



n=l 



H„,H* 



exp [J2 QL + H* n ^n + H n ta* Ln ■ ]Je 



v n=l 



n=l 



It is also useful to write the last L products as a non local action along the boundary 
line 



n 



n=l 



L-l L 



e UnHn = exp 



^ U n H n + U rnU n H m H n . 



(32) 



,n=l 



m=l n=m+l 



4-4- Fermionic action of the PF 

Putting equation (|3*2|) into (|2l)|) . we obtain the Grassmannian representation of the PF 
Q[h] = J Va*VaVb*VbVH*VH expS[a,a*,b,b*,H,H*}, (33) 
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with the action S defined as 

L L L 

S = i^rm + a mnO*mn + ^mnKir) + ^ H n^n + ^ #r 



ta*, 



mn=l n=l n=l 



(34) 



+ 2_j u mU n H m H n + H n H^. 

m<n n=l 

This action can be separated into three terms 



with 



S — 5buik + Smt + iSficid, (35) 

L 

^bulk = 2J (^mn + a mnO* mn + ^mn^r, 



J mn u mn /i 

mn=l 



Afield = u m u n H m H n + ^ H n H*, (36) 

m<n 71=1 

5 in t = 2j ^n L n + H n ta* Ln . 
n=l 71=1 

The PF written as (|33|) is just a Gaussian integral over the set of variables 
(a,a* ,b,b* , H, H*). If we first integrate over the variables (a, a* ,b,b*) corresponding 
to the action <Sbuik + Smt, we obtain a new Gaussian action depending only on fermionic 
magnetic field (H n ,H*). This new action is very similar to that for a one dimensional 
problem. Actually the way we integrate (j33J) is close to solving a ID Ising model with 
a magnetic field and Grassmann variables. In the next section we present briefly this 
case since we will use similar tools later. Our method can then be checked using the 
transfer matrix techniques. 

5. ID Ising model with a homogeneous magnetic field 

The treatment is similar to the 2D Ising model, except that there is only one kind of 
link variables and no mirror symmetry involved. The exact solution in the case of a 
homogeneous magnetic field u = tanh(/3/i) using the transfer matrix method is simply 

2 L Q 1D (h) = (l+t+ v / (l-t) 2 + 4tM 2 ) L + (l + t- ^{l-tf + Atu 2 ^ . (37) 
If we apply the Grassmann transformations as before, we can write an equation similar 

to pa]) : 

Q m {h) = J Va*Va VH*VH exp(<S b uik + <S int + <S fleld ), (38) 



with 



L 



<Sbuik — ^^(a n a* + ta*_ 1 a n ), <S int — H*h n , 

n=l 71 

Afield = 2^ ^nH* + u 2 H m H n , L n = a n + ta* n _ v 



m<n 
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The Grassmann variables a n and H n are anti-periodic in space, and can be Fourier 

transformed using a„ = -4= V r n 1 a„,\ with r a = e 2lnq ^ L . In the new basis, the actions 
& " VL ^i 9+5 1+2 i 

are almost diagonalized 

L/2-1 

Sbuik + Sint = 53 C 1 ~ ^ 9+ i)a 9+ iaJ + i + (1 - tf q+ i)a_ q _ia*_ q _i 

9=0 

L/2-1 L/2-1 

+ E ^ + iK+| + * f ? +i a - g -p+ E H - q -^ a - q ^+tr q+h a* qH ), (39) 

where the bar defines the complex conjugate. For L odd, we have to take care of the 
momenta on the diagonal of the Brillouin zone: it leads to additional terms in 
however, these terms are irrelevant for large L. In the following, we restrict ourselves 
to L even. 

The non local sum in iSfi e id can be written in the Fourier modes as 

J2 H m H n = J2 ^ZZi H q +\ H - q ^ (40) 

m<n q 1+ 2 

L/2 ~ 1 ( 1 1 \ 

= E ( r— i - —731 J E i+\ H -i-\- 

q=0 \ 1+2 1+2 / 

We can separate the previous sums (jHEJ) into independent blocks of four Grassmann 
variables (a q+ i , a* +1 , a_ q _i , a*_ _i). To compute the individual block integrals, we use 
the general Gaussian formula: 



da*da db*db exp(aaa* + abb* + v a + v*a* + wb + w*b*) = aaexp 



ww vv 



a a 



where (v,v*,w,w*) are Grassmann variables and (a, a) two independent complex 
numbers. We obtain 

'L/2-1 



Qm{h) = Qid(0) J VH*VHexp ^ H q+ i_H* q+1 _ + H_ q _i_H*_ 



i-h 



q=0 

+^ H -i-^ H U + ) > ( 41 ) 



where 



2itsin0„,i 

s 1D (o) = n 11 - iffj - r^r^ . («) 

9=0 9+2 



<* 9+ i = -i cot(- Vi), ^9+i = — ( ? + o J • ( 43 ) 



2 q+ i 1+2 l 

Qid(0) is the PF in zero field, and is equal to 1 + t L in this case. 
The remaining integrals over (H, H*) are easy to evaluate, if we use 



J da*da db*db exp(aa* + bb* + aa*b* + (3ab) = l-a/3. (44) 
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After some algebra and simplifications, we finally obtain 

L/2-1 , , 

Qm{h)= J] ( l + t 2 -2tcos^ + i +4tu 2 cos 2 (-^ + i)J . (45) 

Results (J37|) and (|45|) are equivalent when L is even but are written in a different 
way. This has been checked numerically for finite size systems, and analytically in the 
thermodynamic limit. 

Using Grassmann variables to express the PF in terms of products over Fourier modes 
is of course less efficient in the ID case. However in 2D, the method is very similar and 
leads to a final expression which is similar to (J45|) and (|41|) as seen in the next section. 



6. Explicit calculation of the PF 

In this section, we perform the Gaussian integration of the Grassmannian form of Q 
(JHHJ) : We start by reducing the number of Grassmann variables per site, then we 
integrate over the variables in the bulk, in order to obtain a ID action, expressed in 
terms of the fermionic magnetic field. Finally, this last integral is evaluated, leading 
to the PF on the lattice. In order to validate our method, we check that our result is 
identical that obtained by McCoy and Wu, in the thermodynamic limit. 



6.1. Reduction of number of Grassmann variables per site 

In the 2D case, we can easily integrate half of the Grassmann variables, for example 
), by using the identity 

dbda e ab+aL+bL = e LL . (46) 

Since the original measure is do*dfeda*da, moving db to the right of da* implies a minus 
sign. After integrating over (a, 6) we define c mn = a* mn and c* mn = —b mn , which removes 
the minus signs from the new measure dc*dc. We thus obtain the following actions 

L 

^bulk — ^ ^ C mn c m n ^( C mn c mn)( c m-ln ~ C mn-l) ~ ^ C m-ln c m n-li 
m,n=l 
L 

S in t = tH nC0n + {tc[ n _ x + c\ n + C ln )H*, (47) 

71=1 

L 

Afield = 2_j U, m U nHmHn + 2_j HnHn- 
m<n n=l 

Taking into account the different boundary conditions, we Fourier transform these 
variables as in the ID case 

^ L-l L-l 
Cmn = -jr ^ r p r q+ iC pq+ i, c mn = — f p f q+ iC pg+ i, (48) 
p,q=0 P)?=0 
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P 

A 



(0,0) 



o — o 
-o — o 
-o — o- 
-o — o 



(L,L-1) 



q + 1/2 



Figure 2. Integration domain for the Fourier modes in the case where L is even 
(L=6). The set Si of modes corresponds to the white points. The other modes (set 
S2) are obtained from set Si by the symmetry (p, q) — > (L — p, L — q — 1/2) modulo L. 



and compute the bulk action 



L-l 



<S bu i k =^2{l-tf P - tr q+ i - t 2 f p r q+ ^j c M+ ic; + i + 

p,q=0 (49) 
-tr p c pq+ ic L _ pL _ q _i + tr q+ ic* pq+ iC* L _ pLq _i. 

The result implies a coupling between modes (p, q + |) and (L — p, L — q — ^) ~ 
(— p, —? — §)• As in ID, the sum can be expressed as a sum over independent blocks 
containing the variables c_„ , 1 , c _ _ 1 , c* 1 and c* 1 . These different blocks are 
indeed independent if we restrict to momenta (p, q) corresponding to the white points 
(set S\) of the Brillouin zone in figure El In this case, the modes (p, q + ~) belonging 
to 5*i and the modes (L — p, L — q — |) (corresponding to S2) completely fill a Brillouin 
zone. This avoids counting the Grassmann variables twice. We also need to make the 
action more symmetric, which can be written as 

( a pq+h C Pi+k C p<i+^ + ®-p<i+\ c -p-<i-\ c -P-q-\ + fip c pq+\ c -p-q-^ 



r 1+2 PQ+^ -P-I-^J 



(50) 



a pq+ i = l-tr p - tr q+ i - t 2 r p r q+ i, (51) 
P P = t(r p -r p ), (52) 
and for the interacting part 

P,qeSi (53) 

^v-,~^p e %-\ - + c %- q -fp( l + ^+§KV 

The last action, <Sfl e id, is the same as in ID. 
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6.2. Integration over bulk variables : ID action 

The integration over the variables c and c* is performed using the following identity, 
similar to formula (|41j) : 

da*da db*db exp(aaa* + abb* + (3ab + (3a*b* + av + bw + a*w* + &*«;*) 
1 



(«« — exp 



aa — j3(3 



(av*v + aw'w + (3vw + /3v*w*) 



(54) 



We obtain 

£>c*£>c e ,Sbulk+ ' Sin 



n 

P,<?eSi 



- exp (t m+ |^: ? _i^; + i + K q+ \H q+ , j 



where we have identified the different coefficients of the Grassmannian fields forming a 
quadratic action: 

1 1 

1: 



pq+ * L a 



i rv i 
M+5 M+5 



, — ai„„ , 1 ( 1 + 1 ) + <X„_i_ 1 1 1 + tr„ , i 

f3 p /3 +1 V P?+2 V 9+2 ; P9+a V 9+2- 



+^(l + tr g+ i)(l+tf ?+ i)+^ + i) 



(56) 



and 



1 



^ w+ 2 L a m ,ia„. , i 



2 



7 ( tr pPq+ 



i +tr p a M+ i(l + *r„ J _i 



q+i 



pq+h a pq+ 5 fipfiq+ \ 



(57) 



Inside the product (JS3j) . the quantities in front of the exponentials can be simplified 
using cosine functions, 

and are invariant under the transformation (p, g + |) — > (L — p, L — q — |). In (|53j) . the 
product over Si of these coefficients define a bulk PF <2 : 

L-l 

Q o= II [(1 + i 2 ) 2 - 2t(l - t 2 ) (cos0 p + cos0, + i)] . (58) 

p,q=0 

In the thermodynamic limit and in zero field, the free energy per site corresponding to 
<2 is equal to the one corresponding to <2[0], since the boundary conditions do not play 
any role on the bulk properties. In this limit, the second order phase transition occurs 
at a temperature given by the solutions of the equation (1 + 1 2 ) 2 — 4i(l — t 2 ) = 0, or 
t c = \pl — 1, when the cosines, in the long-wave length limit (p, q) ~ (0,0), approach 
unity. In this case the free energy is singular. 

The previous coefficients l pq+ i are not symmetrical in (p, q) since the model itself is not 
symmetrical in both directions. However they are antisymmetric ^p^-i-i — — Tl— pL— q — -• 
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This is useful in order to reduce the summation over the variables (p, q) in the action 
()55|) . This implies 

L/2-1 

V 7„„,iF* ,r, = V 7 2 f!i7* ,Jf xl) 



with 

L — !. 

y 2D i 

p=0 

After simplification, we find that 

1 sin 6 q+ 1 

7p<?+ 5 = 1(1 + t 2 ) 2 - 2t(l - t 2 )(cos# p + cos# g+ i) 

and 

L_1 2ztsin# 



1 z«tsinw 9+ i 
' = ~L (I + t 2 ) 2 - 2th - t 2 )(cos9 n + cos 0^0' ^ 



,,2D _ 

p=0 V 7 V /V p <?+; 

The factors 7 2D play the role of Fourier coefficients of an effective interaction between 
the boundary spins in the magnetic field. We can notice from (JIT]) that coefficients 7 1D 
describe the nearest neighbour interaction of the ID Ising model. Here the spins on the 
boundary can be mapped onto a ID model. 

The factors A pg+ i have a different symmetry, X pq+ i = \_ p _ q _i, which allows the same 
kind of manipulation as before. Hence, defining 

L-l 



p=0 



t ^ (1 - t 2 ) cos(# p ) - t(l + 2t cos 6 q+ i + t 2 ) 
Lj^ (l + t 2 ) 2 -2t(l-t 2 )(cos£ p + ^ +| ) 



(60) 



we can write 

L/2-1 

EA„, i H„, i H* , i + A__ , i H „ i H* i = > A_ , i H„, i H* , i + A „ i H „ i //"* 

P,qeSi q=0 

Moreover, further reduction of the terms containing € pq+ i in the action leads to the 
following simplification 

L/2-1 

p,qeSi 9=0 

Finally, the problem of the boundary field is reduced to a ID Gaussian action in 
Grassmann variables which is given by: 

L/2-1 

(61) 
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with 

This action has the same form as that for the ID problem ([39)1. except for the additional 
terms A and — £ 2 7 2D , which are not present in the ID case. This is due to the fact that 
<2 is n °t the true action in the zero field case: If we integrate (JoTj) with respect to the 
Grassmann fields when u = 0, this will lead to a non zero corrective factor in front of Qq, 
which is however irrelevant in the thermodynamic limit (the free energy corresponding 
to this factor is of order L instead of L 2 ). Q Q is therefore not the finite size zero field 
partition function for the periodic/free spin boundary case. The factor comes from the 
free boundary conditions that restore these conditions. 

Using this ID Gaussian action, the partition function for the 2D Ising model with an 
inhomogeneous or random boundary magnetic fields reads : 



Q[h] = Q J diTdtfexp (Sid) . (62) 

In the following, we will compute the remaining Gaussian integrals in two special cases 
where Aqg/f-u] simplifies. The first one is the case of an homogeneous magnetic field, 
where A gg /[w] = u 2 5 q+ i/ 2 S(q + q') with S(q — q') the Kronecker symbol, and S q+ i/ 2 defined 
in ()43j) . The second case corresponds to the simplest case of inhomogeneous magnetic 
field, when half of the boundary spins is subject to +H and the other half to —H. 



6.3. Expression for the partition function and thermodynamic limit 

In this section, we are interested in the special case of an homogeneous BMF. The 
previous action (|6Hj) reduces to: 

q=0 

L-l 

7 2 ° 1 /T 1 H\ 1 -t^^H^iH „ J +u 2 iH a+ iH a i, 

'9+2 -1-2 1+2 1+2 q+ 2 2 J Z-^l 9+2 1+2 ~?~2' 

with S q+ i defined in (JJ3J). The successive integrations over the blocks of Grassmann 
magnetic fields are easy to perform and we obtain: 



L/2-1 

Q(h)=Q n z qH ( U ), 

q=Q 

with 

Z q+h (u) = (l + A^Kl + l,.^^^ -*\+|)" ( 63 ) 
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We factorize the previous expression, in order to distinguish between the boundary effect 
in zero field and the contribution due to the magnetic field: 

L/2-1 

Q(h) = QoU [(l + A ?+ i)(l + A-,-i)-i 2 (7Si) 2 ] 1 + 

q=0 

The total free energy is therefore written as 



V 2 ft l <"V 2D 



(1 + A + ,)(1 + A_ a ,)-^( 7 2D )2 



*■<?+ 



q+i 



F(h) = —Lk B T In cosh(/%) - k B T In Q + Lf h + Lf t 



field j 



(64) 



with /b a corrective free energy that is needed to restore the free boundary conditions 
in the direction transverse to the field, 



L/2-1 



i E - [< 

5=0 



l + Vi)(l+A^)-t 2 (7? + D .) 2 



(65) 



and /geld the magnetic contribution to the free energy 



L/2-1 



Pfi 



field 



L 



E'» 



q=0 



2r 2D 



(66) 



This decomposition is in agreement with McCoy and Wu's results [T2"| ITT)] . Indeed in 
the thermodynamic limit, we can use the following identity (a > b), 



2jt 



dO 



2vr J a + bcos8 ^/( a - b)(a + b) 



to obtain 



2- 



d9 r . 



1 



R(9 



q+i 



2tt (1 + t 2 ) 2 - 2t(l - t 2 )(cos6 p + cos0, + i) 
with the function R defined by 

R{9) = [(l+t 2 ) 2 + 2t(l-t 2 )(l-cosfl)] [(l+t 2 ) 2 -2t(l-t 2 )(l + cos^)] . 
Then the following coefficients can be evaluated in this limit: 



2it sin 8„ 



q+i 



1 i(l + t 2 )(l-2tcos^ +i -t 2 

1 ^ 2 

2 2 



(67) 



'm+^ 

Using the previous results and after some algebra, we recover the result of McCoy and 
Wu, and the continuous form for the boundary free energy depending on the magnetic 
field is 

4w 2 t(l + cosl, , 

• (68) 



/5/fieid = ^ / d01n| 14 



[l + t 2 )(l-2tcos6 -t 2 ) + ^RjO) 

The expression (jfifij) of the free energy allows the numerical computation of the specific 
heat, even for L small. The results are presented on the figure l6~31 for various magnetic 
fields. 
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k B T/J 



Figure 3. Specific heat for various values of homogeneous magnetic field on the 
boundary for L=20. The inset in a zoom in the region of the maximum. The vertical 
dashed line is the position of critical temperature for 2D Ising model in zero field, 
ksTc/J ~ 2.26. See also the same field dependence in reference |23> figure 3. 



7. Boundary two point correlation function and magnetisation 



In this section, we compute the boundary spin-spin correlation functions along the 
boundary line between two different sites (1, k) and (1, 1), k < I, in the case where the 
magnetic field is uniform. This is the easiest case since we can use simple properties of 
the Grassmann magnetic fields. Using ()14|) . we have 

(<7ik<7i,) cx Tr [a lk aue- m ] , (69) 

a 

L L 



Tr 



We then write (m ^ 0): 

oifcfl + ua lk ) 



m,n=l 



u(l 



n=l 



(70) 



and introduce local magnetic fields u n = u + (it 



u){S kn + Sin) so that 



(a lk a u ) oc m 2 Tr 



I [ ■ + u n a n 



m,n=l 



n=l 
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We remarks that this expression is the PF of the 2D Ising model in the particular case of 
an inhomogeneous BMF (JHJ). The integration over the bulk variables is not affected by 
this change. The difference appears only in the non local coupling between the fermionic 
magnetic fields H n , 



m=l 



cxp 



L-l 



71=1 



m=l n=m+l 



Using the expression for u n , we obtain 



^ u m u n H m H n = u 2 ^2 H mH n + (1 - u 2 )(H k L k + H L Li) + 

L fc-1 



i-M 2 



-H k H h 



(71) 



(72) 



n=k 



n=l 



and therefore, 



,uH n 



\ + {l-u 2 ){H k L k + H l L l ) + 



u 



H k H, + (1 — u 2 ) 2 )H k L k HiL 



k ij k ni^i 



n=l n=l 

Then the two point correlation function can be simply expressed with correlation 
functions (H k Hi), (H k L k ), and {H k L k HiL t ) 

(crikCTn) =u 2 + u 2 (l - u 2 ) (H k L k + H1L1) + (1 - u 2 ) (H k H t ) 

+ ^(l-uYiHkLkHiU). 

The correlation functions involving a product of four Grassmann fields can be written 
in term of products (H k Hi) according to Wick's theorem: 

(HkLkHtLt) = (H k L k ) (H1L1) - (H k L{) (H t L k ) - (HM) (L k U) . (73) 

Using a Fourier transformation, the two field correlation functions are expressed, using 
the usual definitions, as 

L/2-1 



1 > sin 



(74) 



g=0 



,2D 



2^ %+, 



L ^— ' Z n , 1 (u) 



sin 



Each term on the right hand side of equation (J73|) can then be evaluated using the 
previous result: 



,2D 



(H k L k ) (HiLi) 
(L k L t ) 



0„ 



L U ^+§h 2 



(75) 



1 



L/2-1 



*'7 +i cos[0,.i(*- 1 + 1/2)] 



'9+ 



L/2-1 



2D 
9+5 



sin[^i(*-0] 
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The magnetisation can be computed the same way. Using the identity (|70|) we obtain 
(au) =u + u{l- u 2 ) (H k L k ) . (76) 
The connected correlation function is then 

(Wii) - Wik) fa) = (1 - u 2 ) (H k H t ) - u 2 (l - u 2 ) 2 {{H k U) (H t L k } + (H k Hi) {L k U)) . 

These correlation functions, and particularly (H k H[), can be extended for the study of 
more complex configurations of the boundary magnetic field. Without field (u = 0), the 
correlation function between two spins U\ k and o\\ are simply the correlation function 
between the two Grassmannian fields H k and Hi and is given by equation (|74j) . We can 
extract from this the dependence of the magnetisation per spin m in the thermodynamic 
limit near the critical point T c . Indeed, it is usual to define m 2 as the limit of the two 
point correlation function for large separation r = \k — l\: 

lim (a lk au) = m 2 . (77) 

r—roo 

To obtain the main contribution from (|73jl near T c in the thermodynamic limit, we use 
the expressions (JoTj) to compute Z q+ ±(0) in (j^jj) and then make an expansion around 
t c in the low temperature limit (t > t c ) of the different quantities. We first write the 
correlation function as an integral 

2t r ^sin(6) JRC9) . ,„ , 

W=- dg v >y \ > sin(gr); (78) 

1 r , n 2 



S(6) = - ^R(9) + (l+t 2 )(l-2tcos9 -t 



+ 4t 4 sin' 9, 



and expand R and S for At = t — t c and 9 small, which is the region where the main 
contribution of the integral comes from. For R, we find the following expansion: 

■ t c (i-tp - 2+ ; 

2(1 + t c + y/2) 



^R{&) = At(l + t 2 c )(l + t c + V2) 



= At(R + R 2 9 2 + ...) 

where we defined 9 At = 9, Rq = (1 + ^)(1 + t c + y/2) and R 2 = t c (l - t A c )/2. For S, we 
obtain: 

S(9) = At 2 (S + S 2 9 2 + ...), 

where S'o and S2 are numerical coefficients evaluated at t c : So = 16(3 — 2y/2) ~ 2.745, 
S 2 = 4tp ~ .118. After some algebra, we obtain the following behaviour near t c of the 
two point correlation function: 

/ \ 2tAt HA XR ° ■ (70^ 

(a^au }~ / dx— — — — — —smx. 79 

7T J S (rAt) 2 + S 2 x 2 

For rAt small enough, the integral is a constant 7iR /2S 2 , and m is then proportional 

to \^At, which gives the mean field exponent (3 = 1/2 for the boundary magnetisation. 

With the presence of a small BMF, at T c , we use equation (|76|) to compute directly the 

magnetisation: 

, V ,rr r . Ut ( 1 + COS 9 ) JRW) , , 

W / X 7 vr 7 ^) + 2t M 2 (l + cos^)v^) V ^ 
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The quantities R and S are easily evaluated at t c : 
R{6) = 16C (1 -cos#)(3 



COS I 



S(6) = 8C (1 



COS I 



6 — COS 0+^2(1 - cosfl)(3 — cos( 



with Co 



17 - 12y/2. When u is small, the main contribution in (jHUjl is when the 
denominator is small, or when 9 is close to zero. In this case we have ^/Rj9j ~ 
and S(9) ~ 8Co$ 2 . Replacing these values in (|8*0|l we obtain in the limit where «< 1: 

1 2t 

-.ulogu. (81) 



ut 

{(7lk) - — 
7T 



d9- 



10 VQ)0 + 2tu 2 ny/Co 
These results agree with McCoy and Wu's paper [T7j and we can extend the method to 
the case of an inhomogeneous BMF, as we will see in the next section. 



8. An example of inhomogeneous boundary magnetic field 

The previous solution in the presence of a uniform BMF (configuration C a ) on one or 
two sides has be presented in different publications (see for example [2], and |33j), 
but the existence of a ID Gaussian action allows us to solve more general configurations 
of the BMF. In this section, we illustrate that considering the simplest case of a non 
uniform field at the boundary: Between sites (1, 1) and (1, 1) we impose a field h, and 
between sites (1, l+l) and (1, L) a field —h (configuration noted C&). This problem could 
be interesting for the study of interfaces [THI E3 EDI EH E3 EH] and to our knowledge it 
has not been solved exactly for arbitrary values of h on a finite lattice. 
Using our method, it appears that the evaluation of the PF Q(h; I) is a simple extension 
of the previous calculation. Indeed, the product (|75|l depending on the field can be 
simply expressed as 



exp 



u m u n H m H n = exp lu 2 ^2 H m H n - 2m 2 ^ H m ^ H n 



\m<n 



m=l n=l+l 



Then we have 



Q(h; I) = Q Tr 



e s ^ 1 - 2u 2 Yl H ™ Hn 



m=l n=l+l 
I L 



= Q(h) 1 - 2U 2 E ( H mHn) S , 

\ 771=1 n=l+l 

since the other terms from the exponential expansion all include the square of linear 
Grassmann sums and are therefore zero. 

Using the Fourier transformation (|74p. we obtain : 



Q(h) t2^ 



L/2-1 i^B 



L t~r 



L ^— ' Z„ , 1 (m) 1 — cos 0„ , 1 

g =0 9+2 V 7 1+2 



(82) 



The free energy difference between the 2 field configurations Cb and C a is positive and 
is equal to 
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- (3f int = In 1 + Au*\ J2 -y^h-, ■ ( 83 ) 

^ L ^ Z q+1 _(u) 1- cos 9 q+ i J 

This term is added to the total free energy and for I = L/2, sin(# + iZ) is simply 
equal to (— l) q . If the transverse direction has an infinite size, we can use the expressions 
(jnZj) to obtain 




2 l^(-D-cot(W^)| (84) 



L U jR{o q+h )z q+h {u) 



(-I)* cot(Vi/2) 



^ (1 + t 2 )(l - 2t cos + i - t 2 ) + 4tw 2 (l + cos0 9+i ) + jR{8 q+ 



In the ID Ising model, for a system of L spins with periodic boundary conditions and 
with a field configuration C& identical to the 2D boundary line with I = L/2, a similar 
result can be obtained, either by the transfer matrix or the Grassmannian methods of 
the section In the first case, we obtain 



«ff = i»|<i-*)' + 8( " Wl 



1 + t - y/(l -t) 2 + Atu 2 ) L + {l+t 
ln[(l-t) 2 + 4tw 2 ] 




and in the thermodynamic limit this leads to 

(1 -t 2 ) 2 



/3/£? = ln -, K \" J o • (86) 

IJmt [(1 -t) 2 + Atu 2 \ V ; 

With the Grassmannian fields, we obtain a different expression but rigorously identical 
to equation (jS5J): 

/ gtu 2 ^ 1 (-l) 9 cot(^ +1 /2) \ 
- R f. 1D = In 1 - v ; v g+ 2 7 ; (o 7 \ 

P/int ^ L ^ l + t 2 -2tcos(Vi)+4tn 2 cos(^ +| /2) 2 j • 1 1 

In the thermodynamic limit and at zero temperature, f^t is equal to 4J for any non 
zero value h, which is the energy difference between the 2 ground states of Cb and C a , 
all spins following the magnetic field direction in both cases. In the 2D case however, 
if we suppose that the size is infinite in the transverse direction, for small values of the 
magnetic field the boundary spins point all in the same direction (imposed by the bulk 
spins, up for example) below a critical value of the field h c = J(l + 4/L) and for the 
Cb configuration. Indeed the interaction between boundary spins and the ones in the 
bulk are strong enough that a small field — h is not sufficient to reverse these spins. 
In this case /i nt = hL and is therefore extensive contrary to the ID case. Above h c , 



23 



the field —h is strong enough to reverse half of the spins that were up and therefore 
/mt = h c L due to the frustrated couplings with the neighbouring spins in the bulk. 
Figure E] represents the free energy f- m t/L as a function of the BMF for different values 
of the temperature. L = 20 and h c = 1.2 J. There is a change of the curve slope 
at h — h c and low temperature curves suggest the previous reversal picture. Figure 
|3] represents the contribution f- m t/L as a function of temperature for various values of 
the magnetic fields. For h > h c the curves saturate at zero temperature to the value 
fint/L = h c as expected. To check the low temperature behaviour of f\ nt , we can 




Figure 4. Boundary free energy fint/L as function of h for L — 20, h c — 1.2, for 
various values of the temperature. Notice the transition at h = h c corresponding to 
the reversing of half the boundary spins. 



perform an expansion of equation ()84j) for h 1 and T < 1. In this case, we assume 
that t ~ 1 and u ~ 1 — 2exp(—2h/k B T). We obtain 

/ i (-l)«cot(0 o+ i/2) \ 

-/3/ int ^ln 1-4*4 E n+V M 2? fl • ( 88 ) 

I L ^ (1 + M 2 ) - (1 -M 2 )C0S^ + 1 I 

The sum inside the logarithm can be computed using equations (|85|l and (fSTf). with 
£ = 1. We obtain 

/ 2(l-u 2 ) L / 2 \ 

-^■ = ln ( (i-^ + (i + ^ )"' 8ftL (89) 

These results allow us also to study quite precisely the effect of an inhomogeneous 
magnetic field on the spins inside the bulk for a finite transverse size system or a fixed 
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Figure 5. Boundary free energy f- m t/L for L = 20, h c = 1.2, as a function of 
tanh(J//csT) for various values of the magnetic field h/J. 



ratio between the sizes of the two directions (square or rectangular system). In this 
case, the previous domain wall that appears for h > h c due to the reversal of half of the 
boundary spins can propagate or jump inside the bulk for sufficiently high temperature 
in order to lower the free energy. This may cause the total magnetisation to cancel if 
this wall separates two regions of equal number of opposite spins. An exhaustive study 
will be published in a forthcoming paper. 

9. Conclusion 

In this article we have presented a generalisation of the Plechko's method to the 2D 
Ising model in the case of a inhomogeneous BMF. We showed that for any configuration 
of the BMF, this model can be mapped onto a ID model with a Gaussian Grassmannian 
action (J61|) similar to the simpler case of the ID Ising in a magnetic field. Results have 
been obtained for the homogeneous case, in order to validate our results with the ones 
obtained by McCoy and Wu[TT)]. and an extension is made to a special inhomogeneous 
case where an interface develops from one side of the boundary. A simple expression of 
the boundary partition function in term of a ID action is given. Contrary to the ID 
case, the interface free energy appears to be extensive, proportional to the system size L. 
A simple argument gives a zero temperature critical field h c to be the field above which 
a domain wall appears on the boundary and may eventually propagate inside the bulk 
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at higher temperature. This generalisation of Plechko's method for the Ising model with 
non uniform BMF appears to be simpler than methods based on transfer matrix theories. 
Further developments are possible, for example in the study of wetting by a defect plane 
j29| . More generally, it appears however that this method does not answer the question 
of the 2D Ising model with an homogeneous magnetic field in the bulk [TT H ITT ) l3T | I38| . 
even for the 3D Ising problem, since the operations using mirrors symmetries generate 
a Grassmannian action which is no longer quadratic. This is due to the fact, in the 
2D case, that linear Grassmann terms proportional to the field do not commute each 
other easily. The mapping of the boundary region onto a ID action could be used to 
study precisely the properties of boundary random magnetic fields. For the 2D Ising 
case, these random fields are found to be marginal [HH] and logarithm corrections to the 
pure case have been checked numerically jlH]. Analytical computations based on the 
ID general action (jfilj) can be performed using possibly random matrix theory (see also 
jUJ in the context of the Conformal Field Theory). 
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